Trapping Abelian anyons in fractional quantum Hall droplets 
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We study the trapping of Abelian anyons (quasiholes and quasiparticles) by a local potential (e.g., 
induced by an AFM tip) in a microscopic model of fractional quantum Hall liquids with long-range 
Coulomb interaction and edge confining potential. We find, in particular, at Laughlin filling fraction 
ly = 1/3, both quasihole and quasiparticle states can emerge as the ground state of the system in 
the presence of the trapping potential. As expected, we find the presence of an Abelian quasihole 
has no effect on the edge spectrum of the quantum liquid, unlike in the non-AbeUan case [Phys. 
Rev. Lett. 97, 256804 (2006)]. Although quasiholes and quasiparticles can emerge generically in 
the system, their stability depends on the strength of the confining potential, the strength and the 
range of the trapping potential. We discuss the relevance of the calculation to the high-accuracy 
generation and control of individual anyons in potential experiments, in particular, in the context 
of topological quantum computing. 



I. INTRODUCTION 

Shortly after the discovery of the fractional quantum 
Hall effect^, Laughlin^ realized that electrons in such 
a system form an incompressible quantum liquid with 
excitations of fractional charge. These exotic particle 
excitations^'^ are dubbed (Abelian) anyons. To inter- 
change two anyons, one obtains a phase factor e*^ for 
the wave function, where is neither an integral multi- 
ple of 27r as required by bosons, nor an odd multiple of 
TT as required by fermions. The presence of these parti- 
cles with fractional statistics is an indication of topologi- 
cal phases i So far, experiments have confirmed the frac- 
tional charge^iSiii^ but the direct observation of the frac- 
tional statistics remain questionable J^ii^ii^ii^ii^'^^ Re- 
cent experiment s— ^— demonstrated the so-called super- 
periods in the conductance oscillations in a fractional 
quantum Hall quasiparticle interferometer, which appear 
to be consistent with fractional statistics. ^^'^^ However, 
some theoretical worksi^^i^ raised subtleties in the inter- 
pretations. 

A second family of anyons is believed to exist in the 
fractional quantum Hall state at u = 5/2. The even- 
denominator state is belived to be a p-wave paired state, 
known as the Moore-Read state or the Pfaffian state, 
which supports half- flux quantum vortex excitations 
Such particle excitations carry e/4 charge each and, when 
interchanged, not just add a phase factor to the wave 
function, but evolve unitarily in its degenerate (or quasi- 
degenerate for finite systems) ground state manifold. 
They are dubbed non- Abelian anyons, which are also 
speculated to exist at u = 12/5. The existence of the 
non- Abelian anyons, although not confirmed by experi- 
ments yet, is of vital importance to topological quantum 
computing 

In theory, the wave functions of quasihole excitations 
can be written explicitly in analytic functions for both 
the Laughlin case and the Moore-Read case. They are 



also exact eigenstates of some special Hamiltonians with 
short-range two-body and three-body interactions, re- 
spectively. Exact diagonalization of finite systems has 
fruitfully revealed some of these quasihole/quasiparticle 
states. In systems with Coulomb interaction, such 
ground state descriptions appear to be sufficient even for 
electrons on a Corbino disk geometry in Abelian cases^ 
as well as for electrons on a disk geometry in the non- 
Abelian case at = 5/2^ In the latter case, the change 
of the edge spectrum in the presence of an odd num- 
ber of non- Abelian anyons at the origin implies the non- 
Abelian statistics of such excitations. In addition, up 
to four non-Abelian quasiholes have been induced and 
oriented tetrahedrally on a sphere, which results in two 
nearly degenerate states with very similar charge density 
profile (presumably a topologically protected qubit).^^ 

To achieve fault-tolerant quantum computing in the 
topological fashion, one needs to be able to create indi- 
vidual, paired, or a small cluster of anyons. One of the 
simplest experimental approaches is probably to use a 
biased AFM (atomic force microscopy) tip to create and 
trap anyons. One may then easily move the anyons local- 
ized at the tip to realize braiding to fulfill computation. 
However, the feasibility of creating anyons at an AFM tip 
has not yet been systematically studied even on the nu- 
merical level. In a earlier work by one of the authors and 
collaborators,^^ it is demonstrated that a short-range re- 
pulsive local potential (as produced by a sharp AFM tip) 
can induce both +e/4 and +e/2 quasiholes, depending on 
the potential strength, in a z^ = 5/2 system. However, 
a mixture of long-range Coulomb interaction and short- 
range three-body interaction is used, and it is not clear 
whether negatively charged quasiparticles can be created 
in a similar fashion. 

In this work, we study the excitation and trapping of 
both quasiholes and quasiparticles with a local poten- 
tial in a microscopic model of fractional quantum Hall 
droplets with both long-range Coulomb interaction and 



2 



realistic edge confining potential. We focus on the Laugh- 
lin primary filling fraction v = 1/3, although the ap- 
proach can be applied to other filling fractions, including 
the intriguing u = 5/2 case,^^'^^ to obtain similar results. 
We find that both positively charged quasiholes and neg- 
atively charged quasiparticles can be excited generically 
by a finite-range tip potential with appropriate sign and 
strength. We confirm that edge spectrum of the system 
is not affected by the presence of a single quasihole, char- 
acteristic of its Abelian nature. Our results suggest it is 
possible to trap individual anyons, as needed in topolog- 
ical quantum computer proposals. We also discuss the 
stability of anyons when the strength of the confining 
potential varies. 

The rest of the paper is organized as follows. In Sec. 
II, we consider the short-range hard-core potential, which 
generates the Laughlin state and the single-quasihole 
state as exact zero-energy ground states. We consider 
long-range Coulomb interaction in Sec. Ill, where we 
apply tip potentials of (5-function, Gaussian, and expo- 
nential forms. We summarize our results and discuss the 
relevance to experiments in the context of topological 
quantum computing in Sec. IV. 




II. HARD CORE INTERACTION 

In this section, we study the two-dimensional elec- 
tron system on a disk at filling fraction u = 1/3 with 
short-range hard-core interaction between electrons. In 
Haldane's pseudopotential language, Vm = (^i,m- The 
Laughlin state^ at the primary filling factor u = 1/3 
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is the exact ground state with zero energy in the sub- 
space with total angular momentum Mtot = Ml = 
3A^(7V-l)/2 for TV electrons in at least Norb = 3N-2 or- 
bitals. In fact, it is the zero-energy ground state with the 
smallest allowed angular momentum; other zero-energy 
states (for larger Norb) are known as edge states. We plot 
the density profile of the Laughlin state for 10 electrons in 
28 orbitals using palette- mapped 3D plot in Fig.flja) and, 
for comparison, along the radial direction in Fig. [T](b). 

As Laughlin pointed out, the state with a single quasi- 
hole at ^ can be written as 

N 

zi • • • = n i^i - 0*1/3(^1 • • • Zn). (2) 

In the disk geometry, ^ can be placed at the origin to 
perserve rotational symmetry. Obviously, this is a zero- 
energy ground state in the Miqh = 3N{N - l)/2 + TV 
momentum subspace for Norb > ^N — 2. In general, 
there can be additional zero-energy states in the same 
momentum subspace, with the wave function being the 
Laughlin state multiplied by a symmetric polynomial of 



FIG. 1: (color online). The density profiles of the 2D elec- 
trons with hard-core interaction. The Fig. (a) and (b) are the 
density of Laughlin state. Fig. (c) and (d) describe the quasi- 
hole state while we add an external potential WqCqCq at m=0 
orbit (Wo = 0.1). Fig. (e) and (f) are for the quasiparticle 
candidate state, which looks in the same shape as in ref.— . 



order N. Such degeneracy can be lifted either by limiting 
Norb = 3N—1^ or by the addition of an impurity potential 
Hw = WcqCq at the m = orbital. In Fig. (TJc) and 
Fig. m^d), we plot the density profile of the quasihole 
wave function. A density deficiency around the origin 
is clearly visible, indicating the presence of a quasihole 
roughly the size of one magnetic length Ib- 

On the other hand, the quasiparticle state of the cor- 
responding Laughlin state is of some ambiguity. There 
is no zero-energy state obtained in the exact diagonaliza- 
tion at Mtot = Miqp = 3N{N - l)/2 - TV, bacause the 
excitation gap of the Laughlin liquid is finite. Here, it 
is hard to compare it with the variational quasiparticle 
wave function 
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proposed by Lauglin, which is not known as the exact so- 
lution of any simple Hamiltonian. We assume that, like 
the quasihole state, the quasiparticle state is the ground 
state of the Hamiltonian of interested at the appropri- 
ate angular momentum Miqh. Here we plot the density 
profile as such a candidate for a quasiparticle state in 
Fig. me) and Fig.^f). 
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FIG. 2: (color online). The accumulated difference of 
the electron occupation number Xll^o ^^(0 between quasi- 
hole/quasiparticle state and the Laughlin state which oscil- 
lates around 1/3 (-1/3) for quasihole (quasipaticle) state, indi- 
cating the emergence of a charge +e/3 (-e/3) quasihole (quasi- 
particle) 



We plot the accumulated difference of the 
electron occupation numbers Xli^o ^^(0 — 
Si^o [^^^'^^(^) ~ ^^(0] between quasi- 

hole/quasiparticle state (with electron occupation 
number n^^ or n^^) and the Laughlin state (with 
electron occupation number n^) in Fig. [2l The dotted 
line in this figure is the average value, i.e., 1/3 (or —1/3) 
for the quasihole (or quasiparticle) state. This confirms 
that there are ±e/3 charged excitations in a Laughlin 
liquid of = 1/3. In the case of a hard-core potential, 
the size of a —e/3 charged quasiparticle (~ 3/^) is larger 
than that of a +e/3 charged quasihole (~ Ib)- 

A similar density-profile plot of the Laughlin quasihole 
and quasiparticle states in spherical geometry has been 
reported in an earlier numerical work.^^ Studies on the 
Laughlin quasiparticle excitations in the disk geometry 
with long-range Coulomb interaction (in the presence of 
neutralizing background charge to be realistic) have been 
absent. 



III. COULOMB INTERACTION 

In this section we study the excitations in a Laughlin 
liquid of N electrons with Coulomb interaction, confined 
by uniform neutralizing background charge (on a disk of 
radius R) at a distance d above the electron layer. We 
use the disk geometry with the symmetric gauge A = 
(— ^), the single-particle wave function (j)m in the 
lowest Laudau level is: 



^m{z) = (2^2^m!)-^/'z^e-l"l'/^ 



(4) 



where z = x -\- iy is the complex coordinate in the elec- 
tron layer. Projected to the lowest Laudau level, the 



Hamiltonian in the second quantization language reads 



mnl 



where c+ (cm) creates (annihilates) an electron at the 
m-th orbital. V^^^ are Coulumb matrix elements 



= j d^ri J c^V20;,+^(ri)0;(r2)^^n+K^2)(?^m(ri), 

(6) 

and Um the background confining potential 



^R^sJ Jp<R \/\r-p\^ + d^ 



(7) 



In order to study the quasiparticle and quasihole ex- 
citations, we include an external local potential Hw, 
created by an AFM tip, for example. So the complete 
Hamiltonian is 



H — He + H\Y. 



(8) 



In the following, we will consider three different forms 
of Hw' (i) a short-range potential at the origin of the 
disk Hw = WqCqCq; (ii) a Gaussian potential Hw = 
Wg "^Zm ^^P(~^^/2'5^)cJ^c^; and (iii) an exponential po- 
tential Hw We Z)m^^P(~^/OcLcm. 



A. Short-range potential at origin 

A short-range potential can be produced by a very 
sharp AFM tip. By sharp we mean the range of the 
tip potential on the 2DEG is smaller than one magnetic 
length, the size of a single-particle wave function in the 
lowest Landau level. In this case, we can model the po- 
tential by Hw = WqCqCo, located at the origin in our 
disk geometry. A previous study^^ has applied the short- 
range potential to create a single +e/4 quasihole and two 
+e/4 quasiholes (or a +e/2 quasihole) in a model of the 
fractional quantum Hall liquid at = 5/2 with Coulomb 
interaction and an edge confining potential. 

To begin with, we apply the same short-range potential 
Hw to the electron liquid at u = 1/3. We present the 
results of a system of TV = 8 electrons in 26 orbitals 
(large enough so that edge excitations have low enough 
energies). The background charge is still confined to a 
disk of R = yj2N jv = a/48, corresponding to the lowest 
24 orbitals, at a distance d = 0.51b above the electron 
layer. We expect the ground state of Laughlin nature 
has a total angular momentum of Ml = 3N{N — l)/2 = 
84, which is found to be right for zero and small Wq. 
When we increase Wq above 0.26 ±0.01, the total angular 
momentum of the global ground state jumps from 84 to 
92, indicating the excitation of a +e/3 quasihole. The 
density profile of the quasihole state is similar to that 
found for the hard-core potential (Fig. [it and d), in which 
the electron density approaches zero at the origin. 
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FIG. 3: (color online), (a) The low- lying energy states (edge 
states are marked by solid red bars) for 8 electrons in 26 or- 
bit als with Coulomb interaction, (b) The energy spectrum 
after exciting a quasihole at the origin by a short range po- 
tential Woclco with Wo — 0.3. 



In Figure [3l we compare the low-energy excitations of 
the system with and without the quasihole excitation. 
We identify the edge excitations, labeled by solid red 
bars, following the approach developed by one of the au- 
thors and his collaborators.^^ We observe that in energy 
relative to the ground state, the edge spectrum looks al- 
most identical with and without the quasihole, implying 
the Abelian nature of the quasihole. This contrast to the 
case of a Moore- Read state, where the presence of a +e/4 
quasihole changes the fermionic edge excitations. The 
number of the edge states (including the ground state) 
are 1, 1, 2, 3, and 5 for AM = 0-4, as expected by the 
chiral boson edge theory.^^^^^ 

The short-range potential is useful to generate a single 
+e/3 quasihole if the edge confinement is not too strong. 
However, since it only affect the local potential at a sin- 
gle orbital, a second quasihole cannot be induced, since 
one cannot deplete more than (on average) 1/3 charge in 
a single orbital in the Laughlin case. For the same rea- 
son, the short-range potential does not support a single 
quasiparticle (charged — e/3), as one can see from Fig. [2] 
a quasiparticle occupies several orbit als, unlike a quasi- 
hole. Therefore, we proceed to study local potentials 



FIG. 4: (color online). The threshold strength Wc for 
the excitation of a quasihole (red dots) or a quasiparti- 
cle (green squares) in u = 1/3 systems with 6-10 elec- 
trons with Coulomb interaction, using a Gaussian potential 



with width s ■ 



with a longer range. 



B. Gaussian-shaped potential 

We now considered the Gaussian-shaped potential 
Hw = Wg exp(— m^/25^)cJ^Cm, or the potential has 
a value Wg exp(— m^/25^) on the m'th orbital. The width 
of the potential is s, while the strength of the potential 
Wg. In the limit of s ^ 0, the Gaussian potential evolves 
into the short-range (5-potential discussed in the previous 
subsection. 

For fixed s = 2 and d = 0.51b ^ we vary Wg to study 
the change of the total angular momentum of the global 
ground state. For example, in a system of N = 10 elec- 
trons in 30 orbit als, we find the total angular momen- 
tum Mtot jumps from Ml 3N{N - l)/2 135 to 



Miqh = 3N{N-l)/2 + N = 145 at Wg = 0.16 ±0.01, in- 
dicating the presence of one quasihole at the origin. For 
an attractive potential, Mtot drops from Ml = 135 to 
Miqh = 3N{N - l)/2 -N = 125 at Wg = -0.13 ± 0.01, 
indicating the emergence of a quasiparticle at the origin. 
Figure m shows the value of Wg at which Mtot changes for 
systems with 6-10 electrons. We find that the threshold 
values for the generation of one quasihole or one quasi- 
particle approaches a constant value of roughly ±0.15 as 
the system size increases. It is not surprising that the 
threshold value is of the same energy scale as the bulk 
energy gap, but it also depends on the detail of the po- 
tential. 

Figure [5] shows the electron density profiles for (a) the 
Laughlin state, (b) a one-quasihole state for Wg = 0.2, 



and (c) a one-quasiparticle state for Wg 



-0.2. The 



density accumulation or depletion at the origin indicates 
the presence of a quasiparticle or quasihole. Compared 
to the case of the hard-core potential, the quasihole is 
slightly larger, while the quasiparticle is smaller with a 
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FIG. 5: (color online). Electron density profiles for (a) 
the Laughlin state for Wg — 0, (b) one-quasihole state 
for Wg = 0.2, and (c) one-quasiparticle state for Wg — 
—0.2. We consider a system of 10 electrons in 30 orbit als 
with Coulomb interaction. A Gaussian potential Hw — 
Wg exp(— m^/2s^)cJ^Cm with s = 2 is applied. 



well-defined peak at the origin. Therefore, in the more 
realistic case with Coulomb interaction (and not too nar- 
row a tip), the quasihole state and the quasiparticle 
state have roughly the same perturbation to the Laugh- 
lin ground state except for the opposite signs, suggesting 
a quasiparticle-quasihole symmetry. 

Since the quasiparticle state we obtain for the Gaus- 
sian potential in systems with Coulomb interaction can- 
not be easily compared with the variational wave func- 
tion [Eq. ([3])], we want to make sure it is not a stripe 
phase, which arises commonly in systems with Coulomb 
interaction. In Fig. [5](c), we find the electron den- 
sity of the quasiparticle state has a large value around 
the origin and edge. A stripe phase oi N = 10 
electrons in 30 orbitals with a somewhat similar den- 
sity distribution and the same total angular momen- 



tum can be represented by a binary string \^sp) = 
1 110000000000111111110000000000) (a Slator determi- 
nant), in which each digit specifies the corresponding 
single-electron orbitals (from to 29) being occupied (1) 
or not (0). Therefore, we wish to answer the question how 
close the ground state with Mtot = 125 is to the stripe 
phase. For this, we plot the lowest four excitation en- 
ergies (energy difference between the lowest four excited 
states and the ground state in the subspace of M=125) 
as a function of Wg in Fig. [61(a). Obviously, there is no 
crossing/anti-crossing between the ground state (which 
we identified as the quasiparticle state) and the first ex- 
cited states as \ Wg\ increases. This is very different from 
the behavior of the next three excited states, which can 
get very close in energy. We further calculated the over- 
laps between the lowest two energy states and the stripe 
state as a function of Wg in Fig. [Hb). While the over- 
lap is increasing for the ground state, it is only about 
5% for Wg ~ Wg = —0.13 when the ground state in the 
M = 125 subspace becomes the global ground state. We 
therefore conclude that the ground state is unlikely the 
stripe state. 

We extend our calculation to a grid on the area de- 
fined by — 2 < Wg < 2 and < s < 3.5 for the Gaussian 
potential VF^ exp(— m^/25^)cJ^Crn. We choose the 
strength of edge confinement to he d = O.bls (stronger 
confinement) and d = 11b (weaker confinement). The 
results for 8 electrons are plotted in Fig. [71 Generically 
we can devided the parameter space into five regions: the 
Laughlin state {Mtot = SN{N — l)/2), the one-quasihole 
state {Mtot = SN{N — l)/2 + A/"), the one-quasiparticle 
state {Mtot = 37V(7V — l)/2 — N)^ beyond one-quasihole 
state (ground states with Mtot > 3N{N - l)/2 + A^), 
and beyond one-quasiparticle state (ground states with 
Mtot < 3N{N - l)/2 - N). We have also done the calcu- 
lation for A" = 10 with similar results, but on a coarser 
grid. In particular, we do observe the ground state with 
Mtot = ^N{N - l)/2 + 2Ar = 155 for A' = 10, con- 
sistent with the angular momentum for a two-quasihole 
state. We do not find any global ground state with 
Mtot = 3N{N -l)/2-2N = 115 (consistent with that 
of a two-quasiparticle state), but at 117. One might 
tempt to speculate this as one of the two quasiparti- 
cles move away from the origin. Nevertheless, in such 
small systems, it is unnecessary and most likely un- 
reliable to emphasize multiple quasiparticle and quasi- 
hole excitations, so we simply mark the regions with 
Mtot > 3A'(A'-l)/2 + A'and AMtot < 3N{N-l)/2-N 
by "beyond IQH" and "beyond IQP", respectively, and 
do not proceed further. 

The main difference between d = 0.51b and d = 1.01b 
occurs along the boundaries of quasiholes, not along the 
quasiparticle boundaries. This, we believe, is due to the 
fact that for a fixed number of electrons the quasihole 
states (not the quasiholes themselves) have larger size 
than the quasiparticle states, thus more susceptible to 
the edge confinement. The difference is more evident at 
smaller s (sharper tips). In particular, a ^-tip can excite 
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FIG. 6: (color online). Evolution of (a) excitation en- 
ergies of the lowest four excited states and (b) overlaps 
of the lowest two states with the stripe phase |^5p) = 
1 110000000000111111110000000000). The system is of AT = 
10 electrons in 30 orbitals with Coulomb interaction (d — 
0.51b) and a Gaussian potential Wg exp(— m^/2s^)cJ^Cm 
with s = 2. 
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FIG. 7: Ground state diagram for the system of 8 electrons in 
24 orbitals with Coulomb interaction, in the presence of the 
Gaussian tip potential with strength Wg and width s. The 
confining charge is located at a distance of (a) d — 0.51b and 

(b) d^llB- 



a quasihole in the case of d = 1.01b ^ but not in the case 
of d = 0.51b for not too large Wg. This, as illustrated in 
Fig.[7J suggests that a finite width s ^ 21b rnay be more 
robust for the excitation of quasiholes and quasiparticles. 



C. Exponential-shaped potential 



In this section we discuss the exponential-shaped 
potential Hw = exp(— m/^)c]^c^. In real 

space, this corresponds to a Gaussian potential V{z) = 
Wge~\^\ /2cr ^ ^Jiicji niay not be too difficult to prepare 
in experiments. Explicitly, by projecting the real-space 
potential into the lowest Landau level, we obtain the ma- 
trix elements 



w: 



,2 



72m+2 



= w: 



CF 



m+1 



(9) 



where the decay length is ^ = l/ln(l -\-l%/cr'^) and the 
effective strength We = VFJ(7^/((J^ + ^b)- Again, 
is the lowest landau Level wave function with angular 
momentum m. 

After applying the exponential potential Hw = 
Z;^exp(-m/Oc]^c^ with ^ l/ln2 (or a = Ib iri 
real space), we are also able to trap a single quasihole or 
a quasiparticle. Again, we consider a system of 10 elec- 
trons in 30 orbitals, with neutralizing confining charge lo- 
cated at a distance of d = 0.51b above the electron layer. 
The electron density profiles for the Laughlin state, the 
one-quasihole state, and the one-quasiparticle state are 
plotted in Fig. [8|). The density profiles look very simi- 
lar to those for the Gaussian potential discussed in the 
previous subsection (Fig. [5]). Like the Gaussian case, the 
quasiparticle state and the quasihole state have roughly 
the same density perturbation (but with opposite signs) 
to the Laughlin ground state. 
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FIG. 8: (color online). The electron density profiles for 
(a) the one-quasiparticle state (VKe=-0.2) and (b) the one- 
quasihole state (VKc=0.28) for the exponential potential 
Hw = Vl^eXlm^^P(~WOcLcm with ^ = 1/ In 2 (or (7 = Ib 
in real space). The system has 10 electrons in 30 orbitals 
with Coulomb interaction. The neutralizing confining charge 
is located at a distance of d = 0.51b above the electron layer. 



IV. CONCLUSION AND DISCUSSION 

To summarize, we study the trapping of quasiholes 
and quasiparticles by a local potential (e.g. induced by 
an AFM tip) in a microscopic model of fractional quan- 
tum Hall liquids with short-range hard-core interaction 
or long-range Coulomb interaction with an edge confining 
potential due to neutralizing charge. We find, in particu- 
lar, at the Laughlin filling faction u = 1/3^ both quasihole 
and quasiparticle states can be energetically favorable for 
the ground state of the Coulomb system for tip poten- 
tials of various shape and strengths. The presence of the 
Abelian quasihole has no effect on the edge spectrum of 
the quantum liquid, unlike in the non- Abelian case when 
fermionic excitations are present. 

Although quasiholes and quasiparticles can emerge 
generically in the system, its stability depends on the 
strength of the confining potential, the strength and the 
range of the tip potential. Experimentally the quantum 



Hall plateau at = 1/3 was found in a high magnetic field 
(~ 15 T).^ In this case the magnetic length Is ^ 70 A. 
Based on our microscopic calculation, we estimate a opti- 
mal range of the tip potential to be 140 A. The size falls in 
the right range of AFM tip size under current technology. 
The Laughlin state in the context of topological quantum 
computing is of less interest due to its Abelian nature, 
although it can be used for topological quantum mem- 
ory. Nevertheless, it is much easier to model in numerical 
studies than the non-Abelian Moore-Read state, and 
the even more complicated Read-Rezayi (parafermion) 
states We expect the results found here can be of help 
for the excitation and trapping of quasiholes or quasipar- 
ticles in the Moore-Read case in future experiments. In 
the Moore-Read case at filling fraction of 5/2, a smaller 
magnetic field ~ 5 T is usually applied. Thus with 
a longer magnetic length we can have even wider tips, 
which should not be a technical challenge. 

With the well-known difficulties of the exact diagonal- 
ization method in highly entangled systems such as the 
fractional quantum Hall liquids, the search for the ground 
states with a few parameters is a time-consuming job. 
The Moore-Read case is even more complicated, since 
the even-denominator state has a smaller excitation gap 
and is competing with stripe phases One might wish 
to develop more efficient numerical methods to approach 
the ground state properties. One development in recent 
years is the application of density-matrix renormalization 
group (DMRG) method— j^j^i^ to the fractional quan- 
tum Hall systems. We implememt the method in the disk 
geometry with results in excellent agreement with ex- 
act diagonalization in small systems. However, we find 
the time to reach convergence (especially near the origin) 
in larger systems is impractically long for the extensive 
search for the ground states discussed in the current pa- 
per. 
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